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ABSTRACT 

The improvement in ultra compact thermal power 
creation is constrained by the unavailability of fluid 
machinery sufficient for supersonic flow 
circumstances. Conservative turbine designs exhibit 
unacceptable performances related to large 
aerodynamic losses and a narrow operation range. 
This paper provides for the first time in the open 
literature the design procedure, and subsequent 
analysis of the turbine performance of a turbine 
satisfactory for supersonic axial pulsating flows, as 
those encountered in pioneering combustors. The 
design comes up to consider the most adverse 
condition, a steady inlet axial Mach number equal to 
3.5. The possible turbine families were classified by 
the velocity triangles and discussed. An elementary 
issue in supersonic passages is to ensure the normal 
shock at the start of the engine is swallowed through 
the turbine passages, namely the turbine passage is 
started. To ensure self starting capability the turning is 
restricted to lower values than in the conventional 
subsonic turbines. The design procedure was based on 
the method of characteristics, converting the inlet 
uniform flow into a vortex flow field, such that the 
adequate deflection is inflicted to the supersonic flow. 
The presentation of the supersonic passage was first 
assessed and then compared to conservative designs. 
The present design procedure and analysis of 
unconventional supersonic turbines provides strategy 
for the plan and optimization of efficient elevated 
supersonic passages, fit to future tightly packed fluid 
machinery. 

INTRODUCTION 

The thermal efficiency of highly developed mutual 
cycles surpasses 62.8%, allowing for a turbine 
opening temperature of 1703 K [1], However, the 


maximum possible efficiency offered by the Carnot 
cycle for the same firing temperature would be 82.8%. 
While further increase in pressure ratios and firing 
temperature gradually enhance the cycle efficiency, 
novel turbine-based thermal plants offer a potential 
leap in efficiency. In particular, the Humphrey and 
Ficket-Jacobs cycles provide a rise of pressure 
through the combustion process. Fig. la and b display 
the entropy-enthalpy and pressure-volume charts of 
the ideal Joule-Bray ton (constant pressure heat 
addition), the Humphrey, and the Ficket-Jacobs 
cycles. The Humphrey cycle is based on constant 
volume heat addition with an isentropic expansion 
and an isobaric heat rejection. In the Fickett-Jacobs 
cycle, the combustion is an explosive process, utilized 
in the 1920s in a Hozwarth turbine [2], or using 
perhaps a special piston-cylinder arrangement, that in 
practice results in supersonic flow conditions at the 
combustor exit. A multitude of research teams have 
suggested the use of pressure gain combustion for 
power plants [3], The use of a constant volume 
combustor, displayed in red in Fig. la and b, allows 
achieving the same time-averaged combustor exit 
temperature at a higher-pressure level (the turbine 
inlet temperature T4 is fixed by the current material 
technology). Conversely, to achieve the same turbine 
inlet conditions of a Joule cycle, the Humphrey cycle 
requires less pressure increase in the compressor. 
Consequently, the Humphrey cycle offers a potential 
surge in specific power and cycle efficiency. 
Nevertheless to achieve the potential gain in such 
pioneering power plant, the turbine efficiency should 
be above a certain threshold dictated by the combustor 
pressure rise and turbine entry temperature. Fet us 
consider that the Joule engine is equipped with a 90% 
efficiency turbine, Fig. lc demonstrates that the 
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required efficiency of the turbine in the Humphrey 
cycle to achieve the same cycle efficiency than the 
Joule cycle, decreases a function of the combustor 
pressure rise. Humphrey cycles operating at large 
turbine entry temperatures, with low compression 
ratio require low turbine efficiencies. For instance, in 
a Humphrey cycle with T4 141500 K, P3/P2145; with a 
combustor pressure rise of 40%, a turbine with an 
isentropic efficiency of 59% would extract the same 
work than a 90% efficiency turbine in a Joule cycle. 

Nomenclature 

> A Area [m2] 

> C Chord [m] 

> E Relative error 

> M Mach number V/ (y.3f.T)l/2 

> R Grid refinement factor 

> 1,2 Inlet and outlet conditions of a stator row, 
respectively 

> 2,3 Inlet and outlet conditions of a rotor row, 
respectively 


a 



b 



c 



Figure-1 indicate the Single closed loop cycles: a) 
Simplified enthalpy enthropy representation; b) 
Pressure-volume layout; c) Turbine efficiency in a 
Humphrey cycle to obtain the same efficiency than a 
Bray ton cycle with a 90% turbine efficiency.(For 
interpretation of the references to color in this figure, 
the reader is referred to the web version of this 
article.) 


Supersonic passage designs 



Figure-2 
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Supersonic turbines have attracted interest from the 
industry since the 1950s due to the high specific 
power they could provide, allowing a reduction in the 
number of low-pressure stages, and thus leading to 
lighter turbines together with lower manufacture and 
operational costs. 

Fig. 2 a depicts three distinct regions in the vane 
nozzle: a converging (subsonic) section, a diverging 
(supersonic) section and a straight section in the rear 
suction surface. To maximize the specific power, the 
degree of reaction is traditionally very low. The 
design of such impulse type rotors can be achieved 
following two different strategies based on the 
method of characteristics 


Above figure shows a passage designed with the 
“comer flow method”. The front suction side is 
curved to cancel the incoming compression waves 
generated along the concave pressure side, followed 
then by a parallel flow. The flow subsequently 
experiences a comer flow expansion, where waves are 
canceled by the concave surface, until uniform 
parallel flow is achieved at the exit. However, the 
design exhibits zero loading in the central part of the 
passage, and the flow turning is limited. 

Possible in take operating regions shown in Figure3 
below, 





State-of-the-art high-pressure-turbine vane passages 
are characterized by a strong area contraction up to 
the throat. Consequently, when conventional vanes 
are exposed to high supersonic inlet Mach numbers 


the diffuser is never self-starting. Standing shock 
waves upstream of the turbine cascade cause large 
irreversibilities and, in particular, abatement of the 
dynamic pressure, 
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An additional problem is the unfavorable pressure 
gradient across supersonic diffusers that may lead to 
boundary layer separation, which additionally limits 
flow turning in the passage. Furthermore, the 
efficiency of a supersonic turbine is seriously 
penalized by the existence of flow separation on the 
suction surface caused by shock-boundary layer 
interactions. Higher loadings, i.e. higher suction 
surface Mach numbers and lower pressure surface 
Mach numbers are desirable for the starting condition 
as well as to maximize the work produced by the 
turbine. However, adverse pressure gradients will 
become much more pronounced within the passage, 
clearly showing an important compromise between 
starting conditions and flow separation. 

The current paper presents an approach to design 
turbine passages with adequate contraction area ratios. 
The performance of these unconventional turbine 
passages was evaluated numerically using the solver 
CFDp p, both in steady and transient operation, to 
study the movement of a normal shock wave across 
the turbine. 

2-D oblique shock wave - 



1 ^- 

The fundamental purpose of the present research is to 
provide guidelines to designers of an unprecedented 
family of turbine passages suitable for supersonic 
applications. In addition, the performance analysis of 
the components is needed to assess the potential 
performance gain of novel thermal cycles operating 
with internal supersonic flows. 


Procedure of Turbine Design 



Above figure 

Displays all the possible turbine configurations 
considering an axial supersonic inlet flow, the 
velocity triangles and the associated entropy enthalpy 
diagrams, Considering adiabatic conditions, without 
cooling, the turbine channel outlet static conditions 
can be evaluated using the compressible flow 
equations, diffusing channel, and a reaction concept 
as sketched in above figure with a rotor nozzle design. 
The work extraction is maximum for a vane outlet 
swirl equals to 671 and rotor turning equals to 1161 
(P2%P3%581), implying no velocity change across the 
rotor passage if there were no shock losses. This 
impulse design (above figure) delivers a high outlet 
swirl, and both stator and rotor passages experience 
diffusion. In the rotor relative frame of reference, due 
to the upstream stator and peripheral velocity of the 
passage, inlet swirl to the rotor is present; hence flow 
acceleration is possible across the rotor passage, 
above figure depicts the velocity triangle of a turbine 
that comprises rotor passages with flow acceleration 
(nozzle type). Such a design would be the most 
benign for the rotor aerodynamics; however the work 
extraction is limited. An alternative solution, with a 
similar work extraction, would be to consider a stator 
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less configuration, as sketched in above figure the 
rotor is diffusing the flow, but such a concept offers 
numerous advantages, particularly a reduction in the 
number of components. In the following the design is 
concerned with either the stator geometry displayed in 
Fig. 5 a and b or the rotor of above figure. 

r 41 } 

i2} 


[Q j stator or rocor passage- by applying rue mass 
conservation principle: 


m y-TfrTtt 
.*1 




W • 1> 'C-V 


111 


m 


Within an axial turbo machine the pitch (g) is constant, 
assuming a constant channel height Nh, the area ratio inlet to 
any station can be expressed according to the following 

expression Eq. (4): 


A a = tf-g-C Wia„\ ^ 4 

w g cascer^) f0S[tf IM , rr i 

Neglecting pressure losses and assuming constant flow proper¬ 
ties,. the mass flow conservation across the passage yields the 
outlet Mach number (Hxjurr .! as a function of the Inlet Mach 

nunibL-j <A^.ni,., )and ihe passage (m-mngas shown by fcq. i dJ r Fora 

stator the equation is terms of the absolute angle (£i IHie „. r )► 

while fora rotor the turning would be expressed in the relative 
faiiie of reference {0**-&*,*): 

{ V - ] .1 \ 


, y_ i i 

( 1-j —^ aukc 1 ^foullrt tOOriauiift IS) 




This value represents the angle through which the 
flow must turn from M% 1 to the required Mach 
number. 


v - ^ '/d r + >+j{ v^r+t)7(rA)*fn -1 





161 


Tlie critical velocity ratio M~ is eKpressed by the following 
relation Eq. '.7); 


M’ = 


( M a ; y -1)/2 \ '* 
\l -f M 2 l f 11/2 J 


17) 


The circular arcs are designed to maintain the vortex 
flow field region, where the flow velocity is inversely 
proportional to the radius: 


V R = constant 


( 8 ) 


Therefore, the pressure and suction side Mach numbers remain 
constant along the vortex region. When normalizing Eq. b: by the 
critical velocity [V rniNV i) and the correspondent radius of the sonic 
velocity streamline (r*) it can be re-written with dimensionless 
parameters (at R-r\ V-V rnrirn i): 


■tf*-t 


m 


Where R* = Rjr* is Hie dimension less radius of the Vortex flow 
field and M r =lis the dimensionless velocity, Once the 
designer specifies ihe values for ir^ and i^ the values of 
M rlfipfr m and are fixed and R iTPP ^ and R u ,, ry / can be 

determined with Eq. (9j. Considering /j N i ^ and flautto JS the in let 
and outlet flow angles respectively, the amount of turning in the 
circular arcs region are defined by Eq. I0) h 


ftkwvf infer " flmki ~ — l 


H hiuff aleM — $autkt ~ I ^eijefer — ) 

Lfuppmnlrt = flitibri ” (Vupjjp.1 ■ VjniffS 

^upptit.nulhff — ftanUfl — I v npprr — V«4rdJn ) (101 
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In order io fully define the circular arcs, (her sutler flow angle 
fiewiFi is determined with Eq. ,5i. for a selected outlet Majcli 
number (M^rirE)- Since the starting and terminal points are known, 
the circular arcs are now fully defined. The design of the lower 
transitional 4 rcs is performed with respect to non-dimensfonal 
astes = xjr' and y* = yfr* r The arcs are designed in a marching 
manner by adding straight wall segments that produce a small 
change in the flow .ingle and the respective var^non in the Much 
number. According to [ IH | the velocity direction f q>) can be written 
in function of the dimensional radius. 


tp= ± const. 


(in 


/(Jf) = v'O’ + f W~ 1 V in " 1 


) +sin "%-1^-r) (12) 


The positive sign in Eq. [U) stands for the expansion lines, 
while the negative sign represents the compression lines. Applying 
a known boundary condition at x' = 0 where fl.s=0 and JT = 
The equation for the expansion lines is defined by Eq. 1 1 i). 


v=^fr\-fm 


(13| 


Table 1 

Inlet conditions to the supersonic passage. 
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Supersonic passages nu tier propernes. 
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Considering now the discretization of the transition arc into k 
segments where 1 each one turns the flow in At 1 degrees it is 
possible tea determine the direction of each one of (hern with 
Eq. (14). 

Vi.i n lr t = v ^- v ln» / n-< lt - l ^ v 0 4 J 

This definition of tp can now be used m Eq.: O); 

/<*!) = 2vi*r-2{i- 1tfv-2v J ^+/tfl£ Wr ) (15} 

Substituting now the two East terms into Eqs. i> anti 12. the 
vatue of becomes only in function of the dimensionle^ 
radius and the incremental How turning is shown in Eq. i U'.' s 

/(ST) = zvm* - 2(fc- i biv- (£v4/-i-n ftr- fi-1) (i si 

Eqs. 12 .md !7 allow the determination the values of If/ and 
Hence for each k increment its coordinates can be determined 
with Eq* (37). 



The sharp airfoils A, B, C, and D, with the throats 
marked with red dotted lines. For in terpretation of the 
references to color in this figure legend, 


ajjttfcf = — (17) 

■^I-ibJrr = nlt , r^ D * , A.iMi‘ff 

Each segment of the transitionai arc is a straight line parallel to 
the vetocicy direction. ip km therefore the slope of the wall ?m. can be 
computed with Eq. 181 

finally the equation of the wall segment can be determined, 
where k varies from 1 to ~(l J rVr)/AV, 

y* = |#*-<*iE™A + I.PiGpr I+fyiwiT!)r»tjnht 09} 

The design methodology of the upper transition a I arc is 
analogous lo the one used for the lower transitional arc. At the 
outlet of the passage (he procedure is reversed to achieve the 
uniform outlet Mach number uVI). The remainder of the airfoil 
design is thereby completed with the use of straight lines, an the 
suction surface, parallel ro chc How direction, 


The above figure shows the geometry of the four 
different passages. In cases A and B, the throat is 
located at the trailing edge, while for airfoils C and D 
the throat is located upstream. The numerical 
simulations demonstrated that for all geometries a 
constant Mach number region exists as defined with 
the design approach based on the Method of 
characteristics. Reducing the turning (from airfoil A 
to airfoil D) we decrease the pressure loading on the 
airfoil and the contraction ratio. In this research the 
airfoil height was kept constant in order to focus on 
the 2D blade-to-blade analysis and ease the practical 
implementation in a stator row (with a given 
diameter) and control the tip radius in a rotating row. 
Consequently, airfoil D was selected to ensure self¬ 
starting of the passage during the testing of the 
cascade from no-flow into the supersonic regime. 
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Trill 3 

I Lubivic passage design chuiacEPftstycs 
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(a) Main flow features in the supersonic backward 
facing step, (b) experimental shadowgraph [19], (c) 
shadowgraph obtained with CFDJ) j) and (d) static 
pressure distribution downstream of the backward 
facing step. 
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Topology of the structural blocks used to mesh a 
single passage 




@ IJTSRD | Available Online @ www.ijtsrd.com | Volume - 2 | Issue - 6 | Sep-Oct 2018 


Page: 413 






































International Journal of Trend in Scientific Research and Development (IJTSRD) ISSN: 2456-6470 



a) b) 




Schlicren images of an expansion fan/shock wave interaction with a Mach 2.5 frcestream 
and the cut-offs oriented parallel to the freestream: a) knife edge* b) tri-colour 
mask, and cj band lattice mask 
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s b 



a) Mach number contours exposed to the inlet 


conditions of Table 1 

b) Distributions of outlet to inlet total pressure ratio. 
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Results and Discussion 

The mathematical Schlieren visualization displayed in 
above figures illustrates the complex shock patterns 
across the turbine passage. The two oblique shock 
waves generated at the leading edge interact with the 
boundary layers of the neigh boring airfoils, and are 
reflected multiple times further downstream. Above 
figures depicts the static pressure distribution along 
the airfoil, which is characterized by the strong 
fluctuations on the pressure and suction side imposed 
by the impact of the direct and reflected shock waves. 

From the above figures shows the distortion in the 
Mach number field generated by the leading-edge 
shocks. Due to the passage area reduction and 
pressure losses, the passage outlet Mach number is 
reduced to 2.15, while the isentropic prediction (Table 
2) stated a value of 3.4. The leading edge shocks and 
their interactions with the airfoil boundary layer cause 
substantial total pressure losses along the turbine. 

Near the end-walls, extremely high total pressure 
losses were registered, in the boundary layers and 
vertical structures occurring in those regions. Moving 
towards the mid-span, total pressure losses are 
reduced, due to the decreased importance of 
secondary flows. P0 outlet/PO inlet reaches its 
maximum mass-weighted average value of 54% at a 
normalized radius (R Rhub)/ (RtipRhub) of 
approximately 0.78. At the outlet section of the 
computational domain, the mass flow-weighted 
averaged value for P0 outlet /P0 inlet is 34%. 

The performance of a stator less turbine was 
numerically assessed considering as a rotor the 
turbine passage previously designed. The relative inlet 
angle, at design conditions, is zero, but we studied 
eight other off-design cases, with a variation in P2 
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781. In all cases we considered identical relative inlet 
total pressure (40 bar), inlet relative Mach number 
3.5, and that the rotor turns at 12,910 RPM, resulting 
in U2%500 m/s. Hence, the total inlet conditions in 
the absolute frame of reference (P01 and T01) 
increase with the relative inlet angle. 

The above figure evaluates the pressure distribution 
on the turbine passage with the aid of Schlieren and 
iso-Mach number contours, at negative incidence (left 
column), design (central column), and positive 
incidence (right column). At negative incidence ((12% 
8) there are two concurrent effects. Firstly, the 
acceleration along the pressure side (in red) is 
enhanced, and consequently negative loading is 
observed in the front part of the airfoil. This is 
coupled with a weak right-running leading edge 
shock. At about 35% of the axial chord the left¬ 
running shock from the neigh boring leading edge 
impacts on the pressure side, increasing the loading 
substantially above the suction side level. At around 
25% of the axial chord the weak right-running shock 
from the other neigh boring leading edge impacts on 
the suction side. By contrast, the strong leading-edge 
left-running shock is reflected on the pressure side, 
and then impacts again on the suction side at about 
65% of the axial chord. Consequently, the pressure 
levels at the suction side approaches the value at the 
pressure side; such low difference creates a weak 
trailing edge right running shock [22], At positive 
incidence angles ((32% 8), the effects at the leading 
edge are reversed, i.e. a steeper acceleration along the 
suction side (in blue) is coupled with a strong right¬ 
running leading edge shock. This intense right¬ 
running shock impacts on the suction side at about 
30% of the axial chord, originating a bubble of 
recirculating flow that is clearly identified by the blue 
region in the iso-Mach contour. This separation region 
is eventually reattached to the suction side, and from 
50% of the axial chord, the flow is again accelerated 
along the suction side. 

Conclusion 

Present investigate on groundbreaking 
thermodynamic cycle's has shown an unacceptable 
performance of the conventional turbine devices. The 
present research has shown that conventional subsonic 
turbine designs are inadequate for supersonic 
conditions, due to the generation of an intense normal 
shock wave at the inlet. The proposed supersonic 
passages were designed using a methodology based 
on the method of characteristics. The turbine passage 


is comprised of three zones which, (i) convert the 
uniform parallel flow at the passage inlet into a vortex 
flow field, (ii) turn the vortex flow, and (iii) reconvert 
it into a uniform parallel flow at the airfoil exit. 
Several airfoil geometries were designed with varying 
exit Mach numbers. 

The planned design method has been assessed using 
three dimensional Navier-Stokes simulations. The 
computational grid was first carefully selected to 
ensure a grid independent solution based on the CGI 
method. The simulations comprised steady and 
unsteady-transient analysis. The results showed the 
ability of the present design to ingest normal shock 
waves, allowing the passage to operate in the 
supersonic regime. Furthermore, the coupled analysis 
of the density gradient contours together with the 
pressure losses revealed the prime source of loss 
attributed to the leading edge shock reflections across 
the turbine passage. The developed design tool allows 
producing turbine passages twice more efficient than 
the current state-of-the-art turbine designs. 

The aerodynamic presentation of a stator less turbine 
was analyzed at several incidence angles. Severe 
pressure abatement is observed due to shock losses 
and secondary flows. At high incidence, a separation 
bubble on the suction side creates significant losses. 
In contrast to subsonic turbines the static pressure 
increases along the turbine. In spite of the limited 
turning, large values of power are extracted. 
Furthermore, the starting phase from stagnation to 
supersonic regime was analyzed with a transient 
simulation, complex shock patters develop during the 
ingestion of normal shock waves, which typically last 
for about 30 ms. Fluid machinery designers are 
usually constrained to operate in the subsonic regime, 
however the present design approach opens avenues 
for the development of revolutionary ultra compact 
power generation concepts. 
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